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Abstract 

We consider the Chalker-Coddington network model for the Integer Quantum Hall Effect, and 
examine the possibility of solving it exactly. In the supersymmetrie path integral framework, we 
introduce a truncation procedure, leading to a series of well-defined two-dimensional loop models, 
with two loop flavours. In the phase diagram of the first-order truncated model, we identify four 
integrable branches related to the dilute Birman-Wenzl-Murakami braid-monoid algebra, and 
parameterised by the loop fugacity n. In the continuum limit, two of these branches (1,2) 
are described by a pair of decoupled copies of a Coulomb-Gas theory, whereas the other two 
branches (3,4) couple the two loop flavours, and relate to an 8X1(2)^ x SU(2)r/SU(2)2r Wess- 
Zumino-Witten (WZW) coset model for the particular values n = —2 cos[7r/ (r -I- 2)] where r is a 
positive integer. The truncated Chalker-Coddington model is the n ^ Q point of branch 4. By 
numerical diagonalisation, we flnd that its universality class is neither an analytic continuation 
of the WZW coset, nor the universality class of the original Chalker-Coddington model. It 
constitutes rather an integrable, critical approximation to the latter. 



1 Introduction 

The transition between plateaux in the Integer Quantum Hall Effect (IQHE) is a quantum critical 
phenomenon, which was predicted theoretically [H [2] and observed experimentally [3] a few decades 
ago. Although experimentally there is no a priori reason to neglect electron-electron interactions, 
it is usually modelled theoretically by noninteracting particles in two dimensions (2d), in a perpen- 
dicular magnetic field and a random potential. Despite the apparent simplicity of this conceptual 
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setup, it turns out to be very difficult to derive analytically the critical exponents of this transition. 
Important progress was achieved by the introduction of a simple network model which retains the 
salient features of guiding centre motion and quantum tunnelling in the presence of disorder: the 
Chalker-Coddington (CC) model 1^. Extensive numerical studies based on the CC model or other 
approaches have led to good estimates for the critical exponents, notably the correlation-length 
exponent v = 2.37 it 0.02 [6] (a larger value v = 2.593 it 0.006 has also been reported [7]). Also, a 
semi-classical argument [8] yields the prediction i/ = 7/3. 

The CC model is also the starting point for several analytical approaches, like the description 
by a (j-model [9], or a mapping to a one-dimensional (Id) quantum many-body system [101 lllj . 
However, from the point of view of critical lattice models, no exact solution of the CC model has 
been found so far. 

The situation is very different for the spin Quantum Hall Effect (SQHE): the generalisation 
of the CC model to SQHE [12] (which we shall call Sp(2)-CC) maps exactly to classical bond 
percolation, where a large class of exponents are known [13]. This mapping of Sp(2)-CC to classical 
percolation was first observed by Gruzberg et al. [H], who used a supersymmetric (SUSY) spin- 
chain formulation. Later on, it was realised [151 [16] that the SUSY lattice path integral maps 
Sp(2)-CC to a statistical model of lattice paths, which are exactly the hulls of bond-percolation 
clusters. Moreover, a number of SQHE physical observables are expressed in terms of percolation 
correlation functions, and this mapping is valid even at the level of lattice models. 

In this paper, we propose a treatment of the original CC model based on the lattice path 
integral. Since the corresponding statistical model involves paths which may pass through a given 
edge infinitely many times, the number of configurations per unit surface is infinite, and the model is 
not directly tractable by exact-solution methods such as Yang-Baxter integrability and Conformal 
Field Theory (CFT). We therefore introduce a truncation procedure, leading to a series of finite 
statistical models, and focus on the first order of truncation. The arising model is a two-colour loop 
model including vacancies, and with loop fugacity n = 0. 

Integrable multi-colour loop models have been known for a long time [T7]. They were originally 
defined through multi-dimensional height models, but they may as well describe coupled copies of 
classical magnetism models, such as the Potts or 0(n) models, and also the ground state of quantum 
loop models. More specifically, in a two-colour, completely packed {i.e. without vacancies) loop 
model [18[ [T9] , new integrable points were identified through a mapping to a braid-monoid algebra: 
the Birman-Wenzl-Murakami (BWM) algebra [20]. In the present paper, we use a similar approach 
on the loop model arising from our truncation procedure, which is a two-colour loop model including 
vacancies. Generalising to arbitrary loop fugacity n, we obtain four critical branches in the phase 
diagram of this loop model. We then study the critical properties of these branches. 

We find that two of these regimes (denoted 1 and 2) correspond to a pair of decoupled Coulomb- 
Gas (CG) theories, whereas the other two (3 and 4) relate to the SU(2)r x SU(2)r/SU(2)2r Wess- 
Zuniino-Witten coset model, for values n = ±2 cos with r € {1, 2, 3, . . . }. We obtain analytically 
two critical exponents: one of them, Xjnt, corresponds to an elliptic deformation of the integrable 
weights, and the other one, -'^(i,i;adj)) is associated to a perturbation of the weight per monomer. 
The truncated, modified CC model is realised by the n = point of regime 4, but this point is outside 
the validity range for the analytic continuation of the WZW exponents. Our numerical study gives 
the estimate v ~ 1.1 for the correlation-length exponent, and dj ~ 1.71 for the fractal dimension 
of paths. This is clearly incompatible with the IQHE universality class, and hence our integrable 
two-colour loop model is only a crude approximation to IQHE. However, the truncation procedure 
may be carried out to higher orders, possibly yielding more accurate, solvable approximations. 

The plan of the paper is as follows. In Section 2, we recall the definition of the CC model 
and its lattice SUSY path-integral formulation, and explain our truncation procedure, resulting in 
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a two-colour loop model. This truncation is compared in detail with the one used in |101 lllj . In 
Section 3, we use a mapping to a dilute braid-monoid algebra to derive the integrable Boltzmann 
weights of the two-colour loop model, as well as the corresponding Id Hamiltonian. In Section 4, we 
identify the four critical regimes of the integrable model and the corresponding CFTs. Numerical 
and analytical support for the identification of these CFTs is given. In Section 5, we examine in 
more detail regime 4, which contains the truncated, modified CC model at n = 0. We discuss the 
analytic continuation of CFT results, and estimate numerically some critical exponents, including 
the correlation-length exponent v. 

The paper has three appendices. Appendix A contains the details of the mapping to the di- 
lute BWM (dBWM) algebra used in Section 3. In Appendix B, we exhibit a lattice holomorphic 
parafermion V'sC-^) i^i integrable model. In Appendix C, we expose the exact solution of a par- 
ticular point in regime 4, which is mapped to free fermions. This mapping provides a valuable check 
on our results, and also gives a proof that the 0(n = 1) loop model has central charge c = ^. 

2 Truncation of the Chalker-Coddington model 
2.1 The Chalker-Coddington model 

The Chalker-Coddington model ^ is a simple lattice model for the IQHE. The latter consists 
of a two-dimensional gas of non-interacting electrons in a disordered medium, subject to a strong 
transverse magnetic field. In the presence of the random potential, the Landau levels are broadened, 
and eigenenergies are of the form E = {k + ^ h^c + Vq, where k is an integer, Uc is the cyclotron 
energy of the electron in the magnetic field, and Vq is a random part. Let us recall briefly the main 
ingredients of the CC model. 

We consider an electron in the eigenstate of energy E. The spatial trajectories of the electron 
over finite time steps At are modelled by paths on the directed square lattice C (see Fig. [1]), and 
the time-evolution operator over At is denoted U. The operator U reads 

U= Ue (g) ^7, , (2.1) 

edge e vertex v 

with two types of factors: 

• On each directed edge e, the operator Ue takes the particle along e and multiplies the wave- 
function by a random Aharonov-Bohm phase exp{i(pe), where the (pe are independent and 
uniformly distributed on the interval [0, 2tt]. 

• At each vertex v, the operator scatters the particle to one of the outgoing edges. In the 
bases (1, 2) and (3, 4) of Fig. [H is represented by the unitary matrix: 

/ tanh/3 l/cosh/3 \ 

V l/cosh/3 -tanh/3 J ' ^ ' 

The parameter /3 measures the distance to the plateau transition at E = Ec = + ^) The 
critical value is f3c = log(l + \/2), and the corresponding energy perturbation is assumed to behave 
as [U 

{E - E,) a (/3 - . (2.3) 

No exact solution of the CC model is known, in the sense that the critical exponents have 
not been determined analytically. However, very good numerical estimates exist for some of these 
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exponents [SlEllT]. In particular, the correlation- length exponent z^cC; defined by the scaling of the 
correlation length 

^ ^ \E - E.l-"^^ , (2.4) 

has been estimated as [6] 

i/cc - 2.37 ± 0.02 . (2.5) 
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Figure 1: (a) Oriented square lattice C for the Chalker-Coddington model, (b) Labelling of the 
edges adjacent to a vertex of C. 



2.2 Path integral representation 

The problem of solving the CC model amounts to the diagonalisation of a random time-evolution 
operator. We want to perform the average over disorder, in order to turn this into a translationally 
invariant 2d classical model. For this purpose, we use the supersymmetric path integral repre- 
sentation [21]. The following derivation is very analogous to what was done by one of us for the 
SQHE [16], and we use the notations of [l6] throughout this Section. 
The Green's function between two edges ei and 62 is: 



G(e2,ei,z) := {e2\{l - zU)-^\ei) . 



(2.6) 



Here 2; is a parameter which plays the role of the energy in the usual Green's function {E — H)^^: 
roughly speaking z ~ e'^, where E is measured from the filled Landau level. We label cl^gr the 
ends of any edge e, with the convention that it is directed in the sense cr ^ ei, and we introduce 
the complex variables 6i(e),6R(e). The Gaussian measure is defined as: 

j [db] (...):= i y d(Re b) d(Im b) exp(-&*6)(. . . ) , [Bb] := ll[dbL{e)][dbR{e)] . (2.7) 



The Green's function can then be written as a Gaussian integral on the biie), bi?,(e) 

/[D6] bL{e2)bl{ei)expAh 



where the action reads 



G(e2,ei,z) 



A 



(edge) 



A 



(vertex) 



/[D6] expA 



^(edge) _|_ ^(vertex) 



z ^ f)L(e)exp(i 

edge e 



E 



E 

vertex v i ^ j 

V 



'ij(e) , 
b*R{ei)SijbL{ej) 



(2.8) 

(2.9) 
(2.10) 

(2.11) 
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and the notation i -s> j means that i (resp. j) is an incoming (resp. outgoing) edge adjacent to 
V. The next step is to express the denominator in (j2.8p as the inverse of a Gaussian integral over 
Grassmann variables /L,_R(e), /L,_R(e): 

G{e2,ei,z) = I [D6][D/] 6^(62)61(61) exp(A + Ay) , (2.12) 



with the measure 



J [df] (...):= / df df exp(-/7)(. . . ) , [Df] := n[d/i(e)][d/^(e)] , (2.13) 

and Af is the analog of A^, with 6, b* replaced by /, /. 

We denote by an overbar the quenched average over the variables (j)e- A useful formula for this 
computation is 

i/"d*exp(„.-* + „-.-) = f (^^j;. (2.14) 

•^'J m=0 ^ ' 



It easy to see, for instance, that G(e2,ei,2;) = 5{ei,e2)- When studying transport properties, the 
main quantity of interest is \G\'^. We write 

1^(62,61,^)12 = y"[D6][D/] 6^(62)61(61) 6^^+^/ X y"[D6][D/] 6He2)6L(6l) 6^^+^J 

= J [D6i,2][D/i,2] 6^1(62)611(61)612(62)6^2(61) e^"+^/i+^^2+^/2 . (2.15) 

Using ()2.14|) . we get: 

|G(62,61,Z)|2 = I [D6i,2][D/i,2] 6ii(62)61i(6i)6l2(62)6L2(6l) X 

ovT^ r /I (vertex) -(vertex)* .(vertex) .(vertex)*] 
exp I yi^j^ -t- /1^2 "T ^/l "T ^/2 J ^ 

n E V^{[^Li(e)6m(e) + /Li(6)/m(e)] [6^,2(6)6^2(6) + /«2(6)/l2(6)] }"^M-2.16) 



The expression (j2.16p for |Gp can be interpreted graphically as follows. Each term in the 
expansion of the product corresponds to a pair of paths (71,72), where 71 respects the orientation 
of the lattice C (forward path) and 72 follows the reverse orientation (backward path). The two 
paths must use each edge e the same number of times mg. Paths configurations are weighted by 
the elements of the vertex S-mathx, and an additional factor (z*z)'"'=. Note that closed loops have 
a vanishing weight, because the bosonic and fermionic contributions cancel each other. 



2.3 Truncation procedure 

In the form ()2.16p . |GP can be viewed as a two-point correlation function in a classical, two- 
dimensional statistical model for two-colour path configurations. No approximation has been intro- 
duced so far, and thus p.l6p is identical to the value of |Gp in the original CC model. The main 
difficulty in evaluating (j2.16p is that the paths 71 , 72 may go through a given edge an arbitrary 
number of times rn-g, and thus the statistical model has an infinite number of degrees of freedom 
per edge. This type of problem is not usually tractable by exact solution methods, so we need to 
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truncate the statistical model to a finite loop model in order to use these methods. This is very 
analogous to what Nienhuis did for the 0(n) spin model [22] on the hexagonal lattice: in that 
context, the spin model with variables (Sj G M", = 1) was formally mapped to a polygon model 
where edges could be used an arbitrary number of times, but the substitution e"^^*'^-' — >• 1 + JSj • Sj 
in the edge interaction led to a finite loop model, while preserving the 0(n) symmetry of the original 
spin model. 

The truncation we propose consists in keeping only the terms of (I2.16P with rUe G {0, 1}, i.e. 
the configurations where each of the paths 71,72 visits an edge at most once. This preserves 
the boson/fermion supersymmetry, ensuring that closed loops still have a vanishing weight in the 
truncated model. This can be seen as follows. In the original expression (j2.1(jp for the action on 
the edges, we can imagine choosing a different fugacity Ze for each edge (so that it now appears 
inside the summation over e.) This does not affect the supersymmetry of the action. On expanding 
in powers of all the z^, the bosonic contribution to a given edge now enters with a factor (z*Zg)™'=. 
Thus our truncation to mg G {0, 1} amounts to keeping only the terms up to first order in the 
expansion of the partition function in powers of z*Zg, and then setting all the Ze = z again.. Note 
that to this order we have either nothing, or a pair of bosons of different fiavours (1 and 2), or 
a pair of fermions of different flavours, propagating along each edge. The supersymmetry ensures 
that each closed loop is counted with weight 0. At this stage it is simpler to switch to a replica 
formulation rather than using supersymmetry explicitly: we have a model with two flavours of 
boson, such that each edge is either unoccupied, or occupied by each flavour exactly once. Each 
closed loop is counted with a fugacity n, taking then n = 0. The vertices are shown in Fig. [2l 

We briefly comment on how higher order truncations would look in this expansion. For example, 
at 0((z*Zg)^) we would have either 2 pairs of bosons of each flavour, or 1 pair of bosons and 1 pair 
of fermions. (We can never have more than one pair of fermions because the Grassmann variables 
square to zero.) Note that in such a truncation we could give such a configuration a weight different 
from (z*Zg)^ and still preserve the supersymmetry. This points to the existence of an infinite- 
dimensional space of possible super symmetric truncations. However in this paper we consider only 
the simplest. 

We denote by |G(e2, ei, the truncated analog of |G(e2, ei, z)^: \G{e2,ei, z)\^j. is given by 
the same expression as (|2.16p . but with the sum running only over mg = 0, 1. Then \G{e2,ei, z)\^j. 
is interpreted as a two-point function in the loop model defined by the loop vertices of Fig. [2] and 
with loop weight n = 0. 

In the original CC model, the parameter /3 in the S-matrix (j2.2p is staggered. It is useful to 
consider an anisotropic version of this, where it takes the value f3 on the even sublattice of C and 
/3' on the odd sublattice. In this anisotropic CC model, the critical line is [3]: 

sinh/3 sinh/3' = 1. (2.17) 

The Boltzmann weights of the truncated loop model are defined in Fig. [2l For general /3, /?' they 
take the values: 

t ,ui,U2,wi,W2,x =1, a, b, , 6^, —ab (even sublattice), 
t' ,u'i,u'2,w'i,w'2,x' = 1, b' , a', b'^, a'^, —a'b' (odd sublattice). 



(2.18) 



where 



a :=z^cosh , b :=z^tanh^/3 , (o tq\ 

a' := z2 cosh-2 /?' , b' := z" tanh^ /?' . ^^'^^^ 



Note that at the isotropic point these weights are 

1, zV2, zV2, ^V4, ^V4, -^V4. (2.20) 
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Figure 2: Vertices of the loop model arising from the truncation of the CC model. 
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10 


Zc{L,L + 2) 


1.029885 


1.030895 


1.031454 


1.031695 



Table 1: Finite-size estimates of the critical monomer fugacity Zc in the model of Fig. [2l The value 
Zc{L, L + 2) is defined as the solut: 
is the effective thermal exponent. 



Zc{L,L + 2) is defined as the solution of Xt{L, I3c, z) = Xt{L + 2, f3c, z), where Xt{L, /3,z) = ^ log ^ 



2.4 Critical properties 

We now discuss the observables of the model. The most important is the mean square Green's 
function between two edges |G(e2, ei, At z = 1 in the untruncated model this gives the point 
conductance. For a system without any open boundary contacts, this is identically equal to one by 
conservation of probability. It is given by the sum over all pairs of Feynman paths going out and 
back from ei to 62, such that each edge is traversed the same number of times in the forward path 
as in the return path, and weighted by the appropriate S'-matrix elements of the CC model. It has 
been argued |23j that the weights for such 'pictures' are all positive. For z — )• 1~ and on the critical 
line ()2.17p we expect a scaling form 



\G{e2,ei,z) 



~ r 



-2X 



|G|2 p 



r{l-z)^/'^f 



(2.211 



where r = |ei — 62], X^q^2 = 0, and df is the fractal dimension of these pictures (whereby their total 

mass M behaves as r f). The absence of the prefactor r I'^i^ is a consequence of the fact that 
\G{z = l)p = 1 for a closed system. 

In the truncated model, we no longer have probability conservation and so the point z = 1 is 
no longer special. Instead, in analogy with other loop models, we expect to find a different critical 
point, at z = Zc, say, such that the average loop length is finite for z < Zc and diverges for z > Zc- 
The point conductance \G{e2,ei, z)\^j. now corresponds to the weighted sum of a pair of black and 



grey paths connecting ei and 62- On the critical line and as z 



we expect the same scaling 



form as in ()2.2ip with (1 — z) replaced by {zc — z), but not necessarily with same exponents as in 
the full model. In Figure [3] and Table [U we show the numerical determination of Zc using the two 
largest eigenvalues Aq, Ai of the transfer matrix. These eigenvalues define the thermal exponent Xt 
through the CFT form of the free-energy gap: 



Ao 2TrXt 



(2.22) 



More generally, it is possible to consider 'watermelon' exponents X^^^i^ corresponding to ii 
black and £2 grey paths originating from the vicinity of a given edge. The truncation constraint of 
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course implies that these cannot originate on the same edge for i > 1, but we imagine taking the 
scahng hmit where edges a finite distance apart on the lattice are mapped to the same point. These 
operators are well suited for a transfer-matrix-based numerical analysis j24j . 

In particular we see that Xi^i corresponds to X^q^2 in (j2.2ip . Also, since z*z counts the number 
of edges connected to 2 black and 2 grey paths, we have 

df = 2- X2,2 ■ (2.23) 

Using transfer-matrix diagonalisation, we obtain the value 

X2,2 =Xt^ 0.3 . (2.24) 

Finally, we evaluate the correlation-length exponent v associated to a perturbation of the pa- 
rameter (3 away from /3c- For the lowest free-energy gap we expect the scaling form 

(2.25) 

The best data collapse is obtained for the value (see Fig. U]): 

(2.26) 



log^ F 



{13 - 13,) 



2.5 Relation to Hilbert-space truncation 

We close this Section by comparing our approach to earlier studies |10^ [TT] of the IQHE problem 
based on a different truncation procedure. Our method consists in writing the lattice path integral 
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(\1 




{(3 - /3e) X L^/'^ 



Figure 4: Data collapse for the effective thermal exponent Xt{L, /3, Zc), under a perturbation of the 
parameter /?. The value used for this plot is l/i^ = 0.9. 



representation for the mean conductance using the supersymmetry trick, and then truncating the 
infinite sum over the paths, to keep only the self-avoiding paths. This gives us the well-defined loop 
model of Fig. O where we will tune slightly the Boltzmann weights to obtain an integrable point 
(see Section [3|). 

In contrast, in [10^ lllj. one starts from a two-dimensional single-particle Hamiltonian including 
Gaussian hopping coefficients, and computes its supersymmetric path integral. The resulting action 
is then interpreted as the action of a one-dimensional many-body supersymmetric Hamiltonian 
Hmb, given in Eqs. (3-5) of [11]. This Hamiltonian is expressed in terms of the coefficients of a 
superspin matrix. In this model, the Hilbert space for each site is infinite-dimensional (each site can 
be occupied by an arbitrary number of bosons). The idea is to truncate this Hilbert space down to 
dimension D, and follow the behaviour of the energy gap as D increases. The model is not critical 
for finite D, but it becomes critical in the limit D — )■ oo. 

Let us shown how to relate the terms of Hmb in the truncated space of dimension D = 5, to the 
generators which encode the loop model of Fig. [2j We first get rid of the (— 1)-' factor in Hmb [TT] . 
This is done through the change 

Cf-j —7- —c-^j for j = 2 mod 4 or J = 3 mod 4 , 

without affecting the (anti-)commutation relations for the bj,Cj. We obtain the Hamiltonian: 



-f^MB = ^ 
where the signs ga are given by 



16 



a=l 



9a 



1 if a = 1,2,10,12,14,16 
-1 if a = 3, ...,9,11,13,15. 



(2.27) 



(2.28) 
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We decompose .ffMB as a sum of generators 



Hmb = 5^ -("),-(J,+e,+/, + (l + r,)[()) ), + ( (0,+2())(0, 
j=i { 



(2.29) 



where we have defined 

((nl)j '■- 



q3q3 i c? c? , clS clS c-16 cl6 



q4 q4 _i_ cS cS _i_ q13 q13 014 014 



(» ), + ())((). 



( (0. + ())((). := Ushl + Sl^ + S^^, + Sl 



^1 



ol 



q5 

-^i+l 



I olO olO I Q12 q12 



Q2 _ iWcS _ 1 
2 ) y^j+l 2 



-S*! + 2 ) ( 'S'i+1 + 2 



q9 q9 



oil oil 



In terms of the creation/annihilation operators, the above generators read: 



(2.30) 



, in) )j 



0) ). + ())((). 
( (0, + ())((), 



(2.31) 



In the D = 5 truncated space, each site is either empty or occupied by two particles of opposite 
spins (t^i)- If each spin is interpreted as a loop color, the above generators (when restricted to the 
D = 5 space) obey a dilute two-color Temperley-Lieb algebra with loop weight n = 0. Hence, they 
represent the vertices ni,n2,wi,x of the loop model defined in Section \T3[ In particular, the ej 
and fj form two decoupled Temperley-Lieb algebras. 

Note that, in this context, the generator for the W2 vertex, Ej = ejfj, cannot be realised by 
a linear combination of the S^Sj^^, but it may be a linear combination of the (SjSj^i)'^. So 
introducing Ej terms in the Hamiltonian leads to higher-order terms in Hmb, and most probably it 
breaks the invariance with respect to the supersymmetric charges Qi,2- However, we have shown that 
the supersymmetric model -ffMBj when restricted to the D = 5 space, corresponds to a particular 
manifold in the phase diagram of the two-color loop model. 



3 Construction of an integrable critical loop model 

In the preceding section, we truncated the the Chalker-Coddington network model to yield a two- 
color loop model that is simpler to analyze. To make further progress, we modify this model 
further. We augment it by allowing the "straight-line" vertices with weight v illustrated in Fig. O 
We also generalize it by allowing the weight per loop n to not only be zero, but to vary in the range 
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n € [—2, 2]. By utilizing the results of |25j . we will show in this section that this modified model 
for all values of n in this range has an integrable line, and includes several critical points. The 
remainder of the paper will be devoted to the study of the critical behaviour. 

When the straight-line vertices are allowed, the loop model can no longer be related directly to 
electron trajectories in a potential. In the original CC model, the 'checkerboard' structure of the 
lattice (or, equivalently, the alternation of arrows on the edges of C) is essential to the interpretation 
of the paths as electron trajectories along the contour lines of the random potential. However, 
several arguments indicate that the truncated but unmodified loop model of Fig. [2] is in the same 
universality class as that of the modified model. In other words, one can obtain the unmodified 
model by perturbing the critical line with irrelevant operators. 



Ul 



U2 









1 











W2 



Figure 5: Vertices of the augmented dilute two-color loop model. 



One argument for the equivalence of the two stems from the relation of this two-color loop model 
to that studied in [19] . There the completely packed version was studied; in the notation used here 
this corresponds to setting the Boltzmann weights t = ui = U2 = v = 0. It was shown that at 
least for weight per loop n > the model has a critical point when x/w2 is tuned appropriately. 
Moreover, at this critical point, numerical evidence strongly suggests that dilution (i.e. non-zero t, 
Ul and U2) is irrelevant. We will provide additional evidence by finding that for certain discrete 
values of n > -v/2, the critical point of the completely packed model and that of the modified model 
studied here are described by the same conformal field theory. Neither of these arguments applies 
when n = 0, but all the critical exponents we have computed {Xt,X2^2,'^) for both the truncated 
CC model and the integrable model at n = in regime 4 (see Table [2]) agree, up to our numerical 
precision. This strongly indicates that the integrable model at n = in regime 4 is in the universality 
class of the truncated CC model. 

In this section, we give the Boltzmann weights of the integrable critical line in the loop model of 
Fig.O These weights are expressed in terms of the generators of the dilute Birman-Wenzl- Murakami 
(dBWM) algebra, so that the solution of the Yang-Baxter equation found in |25j can be used. In 
Appendix A, we review the BWM algebra and its graphical presentation. The braid group can be 
represented in terms of the BWM generators, and can then be used to find invariants of knots and 
links generalizing the Jones polynomial |20j . 

An alternate way of obtaining the Boltzmann weights of the integrable critical line is to search 
for holomorphic observables on the lattice. These are operators whose expectation values satisfy 
the lattice analog of the Cauchy-Riemann equations. This method is described in Appendix B, and 
yields the same weights as those found in [25] using the dBWM algebra. 

3.1 Critical completely packed loop models 

We first review the critical completely packed loop model, arising for example in the Fortuin- 
Kasteleyn expansion of the Potts model [26]. Each vertex of this model has the two possible 
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configurations displayed in Fig. [6l The partition function is conveniently written in terms of the 




Figure 6: Action of 1 (left) and ej (right) on a pair of strands at positions j and j + 1- The transfer 
matrix direction is upwards. 

generators of the Temperley-Lieb (TL) algebra [27]. This algebra for a system of width L has L 
generators ej acting at positions j = 1, 2, . . . , L as well as the identity 1, which obey the relations 

Cj = n Cj , ejej±iej = ej , CiCj = ejei for \i — j\ > 1. (3.1) 

The first of the relations encodes the fact that the weight for a closed loop is n, while the second 
encodes the fact that the weight does not depend on the length or the shape of the loop. 
The Boltzmann weights of the integrable critical loop model are then 

Rj{u) = sin(20 — u)l — sin it , (3-2) 

where n = —2 cos 20 and |n| < 2. The transfer matrix for an even number of sites L is then 

T = -R1-R3 . . . RL~iR2Ri ■ ■ ■ Rl ■ (3-3) 

It is straightforward to use the TL algebra to verify that these Boltzmann weights satisfy the 
Yang-Baxter equation 

Rj{u)Rj+i{u + v)Rj{v) = Rj+i{v)Rj{u + v)Rj+i{u) (3.4) 
and the inversion relation 

Rj{u)Rj{-u) = sm{2e - u) sm{2e + u) 1 . (3.5) 
Braid group generators bj and bj^ are found by taking u — t- ±ioo: 

Rj (ioo) oc bj = e"*^ 1 + e'^ ej , Rj {-ioo) (x bj^ = e*^ 1 + e"'^ ej . 

These satisfy the braid-group relations (1A.1[) and (IA.2j) as a consequence of the Yang-Baxter equation 
and the inversion relation respectively. 

The critical completely packed loop model on the square lattice is in the same universality class 
as what is usually known as the 0(n) model in its dense phase. Well-established results on the 
dense 0(n) model |22] give the central charge of the CFT describing the scaling limit to be 

1 3(7r - 2ef 

'^O(n) = 1 ■ (3.6) 

The Boltzmann weights of the completely packed doubled loop model studied in [191 [18] can 
be written in terms of the generators ej and fi of two independent TL algebras. This model is 
displayed in Fig. [5] with t = ui = U2 = v = 0. In this picture, the ej acts on black loops while the 
fj act on grey loops, while the transfer matrix goes to the northeast. Thus the vertex with weight 
wi corresponds to the generator 1, the vertex with weight IU2 corresponds to Cjfj, while those with 
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weight X are Cj and fj. Since the Cj's and the /j's commute, we have immediately that the Bj,B- ^ 
defined by 

BT^ := (e^n + e-^^ e,) (e^n + e-^Vi) ^ ' 

also generate a braid group. Similarly, TL generators with loop weight N = n? may be constructed 
as 

E,:=e,f,. (3.8) 

Using the relations (j3.ip for the e^'s and /j 's, it is straightforward to show that the Bj, BJ'^ , Ej 
generate the BWM algebra described in Appendix A with parameters N = 'n? = (—2 cos 20)^, 
w = e*®^ |19j . The doubled lines here correspond to the single lines displayed in Appendix A, as 
is apparent by comparing Figs. [5]and[T31 Writing the Boltzmann weights in terms of this algebra 
is useful because solutions of the Yang-Baxter equation involving the BWM generators have long 
been known [28] . Prom this solution, a critical point for the coupled completely packed loop models 
for n > \pl was found |19t [TS] . With the parameterisation 

vr 

n = 2 cos ■ 



r + 2 ' 

in the isotropic case wi = W2, the critical point is at x/wi = Ac, where 

"7r(r - 2) 



Ac = -\/2 



sm 



4(r + 2) 



(3.9) 



At integer values r = 2,3,4..., this critical point was identified with a particular conformal field 
theory, the WZW coset model SU(2)r x SU(2)r./SU(2)2r. This conformal field theory has central 
charge 

3^2 

Cr = 7 TT r . (3.10) 

(r + l)(r + 2) ^ ^ 

For 1/Ac < x/i^i < Ac, the doubled loop model has a critical phase corresponding to two decoupled 
completely packed loop models. The central charge is thus twice (j3.6p . 



3.2 The integrable critical line 

We now can use the results of Grimm and Warnaar [25] to find an integrable model involving all 
the vertices in Fig. [5l We are interested mainly in the critical points, which can be written in 
terms of the dilute BWM algebra. The dilute BWM algebra extends the BWM algebra described 
in Appendix A to include edges of the lattice uncovered by strands. In the two-color loop model, 
these amount to allowing vertices to be empty of both colors. The dilute generators act identically 
on the two colors, and so include the remaining vertices in Fig. [5l In an obvious notation, we then 
can write the ^-matrix as 

R,{^) = %)( ),+ni(v.)[()) ), + ( m+U2{^)[r), + {Jj] 

+ (X)i] + + W2{ip)E, + xi^)Xj . (3.11) 

In terms of the TL generators introduced in the previous section, Ej = ejfj and Xj = ej + fj, while 
Ij takes value on the dilute configurations and 1 otherwise. 

Since the non-dilute vertices satisfy the BWM algebra, it is simple to show that the operators 

Bj:E„I„{ )„()) )„( (()„(^)„(J„(/')„(X), 
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constructed from the two-color loop model satisfy a dilute BWM algebra. Namely, with doubled 
lines here corresponding to single lines in Appendix A, and the Bj defined in ()3.7p . these operators 
generate the dilute BWM algebra with parameters (A^ = {q + q~^)'^ ,u! = q^), where q = e^*^. 

In [25] . an integrable model based on the dBWM algebra was derived. With n = —2 cos 20 as 
before, its Boltzmann weights are given by 



t{^) = 


— cos(2ip — 39) — cos 


59 + cos 39 + cos 


Ul{ip) = 


-2sin26'sin(¥P - 361) 




U2{^p) = 


2 sin 29 sin 




v{(p) = 


—2 sin ip sm{ip — 39) 




wi{ip) = 


2 sin(93 — 29) sm{(p 


- 39) 


W2{(p) = 


2sm(ps'm{ip — 9) 




x{(p) = 


2 sin (p sm{ip — 39) . 





(3.12) 



We denote by (po the isotropic value which is closest to zero: 

The universal properties are independent of the anisotropy parameter ip (as long as (p lies between 
and (/?o)) but depend very strongly on 0, as we shall see. At the isotropic point ip = 39/2, the 
weights can be rescaled to 

t = 2 cos 36^ + 2 cos 29 + 1 

Ui = U2 = 4 COS I COS 9 

v = 2cos9 + l (3.14) 

Wl = W2 = 1 

x = -(2 COS 61 + 1). 

The integrable model defined by (|3.1ip - (|3.12p obeys the following properties: 

• The isotropic weights are invariant under the transformations 9 — t- 271 + 9, and 9 — )• —9, so 
the range of inequivalent couplings is G [0, vr]. Each value of n G [—2,2] appears twice in 
this interval. 

• Since there are no loop ends, the number of loops mod 2 is the same as the number of x 
vertices mod 2. This allows us to change the sign of n by absorbing the sign in the weight 
x: {n, x) —7- (— n, —x). Thus there are four distinct critical points for each value of n G (0, 2), 
while there are two for n = and n = 2. 

• The weights satisfy the inversion relation 

R{ip)R{-ip) = 4sin(26l - ip) sin(26' + ip) sin{39 - ip) s\xi{39 + ip) 1 . (3.15) 



• Rotating by 90° is equivalent to sending ip ^ 39 — ip. 

• The weights are trivial when ti = 0: -R(O) = 2 sin 20 sin 30 1 . 

• The eigenvalues of the transfer matrix are preserved under (ui, U2) — )• (— ui, —U2) and v — )• —v. 



14 



3.3 The quantum Hamiltonian 



To gain intuition into this doubled loop model, it is useful to find the equivalent Id quantum 
Hamiltonian by taking the very anisotropic limit 99 —t- 0. The Hamiltonian is found from the 
transfer matrix TL{ip) for L sites by 



H := 2 sin 26* sin 30 



dlogTLiif) 



dip 



+ 2L sin 56 1 , 



ip=0 



yielding 
H 



^ |4cos 46* sin 6* 



),+2cos20sin30 [()) ), + ( ((),] + 2 sin 20 [(^),- + (J,] 



+2 sin 30 [(y )j + (X) j] - 2 sin Ej - 2 sin 30 Xj | . 



(3.16) 



To find the Fermi velocity Vf, we assume that in the scaling limit this Hamiltonian is that of 
a conformal field theory. In the next Section, we will present much evidence in support of this 
assumption. In a conformal field theory, the ground-state energy (the lowest eigenvalue of H) is 



El ~ Lec 



vrc 
6L 



(3.17) 



where c is the central charge. Let Aj^((/9) be the dominant eigenvalue of the transfer matrix. The 
analysis of Appendix B indicates that the free energy of the loop model on a rhombic lattice with 
angle a is given by [— log A^^ ((/?)] , where a = 7r(/?/(2(^o) and (^0 is the isotropic value, as defined in 
()3.13|) . In a conformal field theory, one expects [29] 

(3.18) 



- log A?^ ((^) ~ L/00 (a) - ^ sin a . 

bL 

Differentiating (|3.18|) around = and comparing with (|3.17|) yields 

27r sin 20 sin 30 



2v3o 



39 
39 



if < < f 

f-vr iff<0<^. 



(3.19) 



4 Identifying the critical theories 

In this Section, we present what we believe is convincing evidence that the doubled loop model 
with Boltzmann weights (I3.12P is critical. We find the presumably exact central charge of the 
conformal field theories describing the scaling limit, and also give some of the dimensions of fields. 
We do this by a combination of calculations exploiting the integrability, comparison to a similar 
integrable model, and exact diagonalization of the transfer matrix and the Hamiltonian for widths 
up to L = 14 sites. 



4.1 The four regimes 

This critical line is parametrized by the value of € [0, vr], related the weight per loop by n = 
—2 cos 20. Since the Fermi velocity vanishes at = ^, ^, and has a discontinuity at = it is 
natural to expect that the physics is discontinuous if is varied across these values. We thus divide 
the critical line into four regimes, as described in Table [2j 
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regime 


t^-range 


parameterisation 


j_ 1 1 

central charge 


1 


< < f 


n = -2 cos ^ 


^ _ o [-1 6 1 1 1 
^ ~ ^ (r+l)(r+2)J ^ 2 


2 




n = 2 cos ^ 


„ _ o 1 6 

L (r+l)(r+2)J 


3 




n = 2 cos 


^ ~ (r+l)(r+2) ^ 2 


4 


f < < vr 


n = -2 cos ^ 


3^2 

~ (r+l)(r+2) 



Table 2: The four regimes of the integrable loop model. 



All known integrable models with Boltzmann weights parameterized by trigonometric functions 
of the anisotropy parameter ip are critical, and this is no exception. One argument for this is the 
existence of the lattice holomorphic operator described in Appendix B. Another is the inversion- 
relation calculation done below, which shows that with standard assumptions about holomorphicity 
in the free energy is singular as this critical point. A numerical check is to use exact diagonal- 
ization to find the largest eigenvalue of T and/or the ground-state energy of -ff, and then fit the 
results to (j3.17p or (|3.18p . To extract the central charge c, we use two different-length systems to 
get rid of the extensive piece Leoo- Doing this, we find the results given in Fig. [71 We see a very 
nice convergence to the critical behavior as expected. 

We combine these results with other arguments to conjecture exact formulae for the central 
charge for all Q. We can also identify precisely which conformal field theories describe some critical 
lines. There are two types of conformal field theories known to describe doubled loop models, and 
both occur along this critical line. Unfortunately, the value of n = at = ^ of interest for 
the truncated CC model lies in one of the regions where we do not understand the conformal field 
theory. As is apparent from Fig. [71 we do know that c = as required there. 

At several special values of 0, the model simplifies. Namely, when n = ±1, all loop configurations 
receive the same weight (if n = —1, we transform (n, x) — ?• (— n, —x) as explained in Sec. 13. 2p . Thus 
when computing the partition function, we can sum up the four completely packed vertices to give 
a single one with weight w\ + 2nx. 

For 6 = ^ at the isotropic point, x = v = 0, so this reduces to a six-vertex model with no 
staggering. Here the usual parameter [30j has value 

a2 + 62 _ c2 



r.2 

A= " ' " =-1. 
2ab 



so this is in the same universality class as the antiferromagnetic Heisenberg model. Thus the central 
charge is c = 1 and the first thermal exponent is = i, in agreement with the numerical results 
in Figs. [THS 

At 6 = ^ and 9 = we obtain a staggered version of the eight-vertex model. Ordinarily 
the staggered eight- vertex model is not solvable, but as we detail in Appendix C, this one is not 
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e 



Figure 7: Numerical estimates for the central charge in the four critical regimes. Different symbols 
represent data points for consecutive system sizes: L = 4, 6 = 6,8 (x),L = 8,10 = 

10, 12 {u),L = 12, 14 (■). Full lines represent the predicted exact values from Table [2j 



only solvable, but can be mapped onto a free-fermion theory. There we show that there are two 
Majorana fermions present, but only one of the two is critical. Thus the central charge is c = ^ 
here, again consistent with the numerics. 

4.2 Computation of an exact scaling dimension 

Since the model is integrable, it is possible to derive some quantities exactly. Here we extract 
the dimension of an operator in the critical theory as a function of 9. This is possible because at 
certain discrete values of 0, there exists a deformation away from the critical point preserving the 
integrability [22]. The inversion-relation method [3D] yields the free energy along this deformation, 
and by analyzing its expansion around the critical point, we extract the value of the exponent 
fint. This then yields the dimension of the operator which when added to the action causes the 
deformation. 

It is convenient to parameterize the loop weight n within each of the four regimes by a param- 
eter r, 

71" 

n = 2ecos , (4.1) 

r + 2 

where e = — 1 in regimes 1 and 4, and e = 1 in regimes 2 and 3 (see Table [2]). 

The integrable deformations resulting in unitary field theories occur at integer values of r in 
all four regimes. Here the dilute BWM algebra admits a "height" or "Restricted Solid-On-Solid" 
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(RSOS) realization (25]. Instead of treating the loops as the degrees of freedom, on the dual lattice 
one places height variables, which are integers restricted to a certain interval. The loops then play 
the role of domain walls separating regions of different heights. 

The inversion-relation method is a way of computing the free energy exactly after making as- 
sumptions about its holomorphicity properties as a function of ip. The free energy satisfies con- 
straints following from the inversion relation (14.20 below, and the fact that sending — )• 30 — 
rotates the lattice by 90°. The holomorphicity assumptions then give a unique solution to these 
constraints. Parameterizing the deformation in our case by p, the inversion relation becomes [25] 

R{p,p)R{-ip,p) = (4p)"i ei{2e - ^,p) ei{29 + v?,^) 9i{3e - ^,p) ei{w + ^,p)i, (4.2) 

where 9i{u,p) is the standard elliptic theta function. This indeed reduces to (j3.15p in the critical 
limit p — 7- 0. From this, it is simple to show that the inverse of the transfer matrix in the diagonal 
direction is given by forming a transfer matrix out of products of Rj{—ip,p). 

Conveniently, both ()4.2p and the behavior under rotational symmetry are identical to that of the 
model studied in [31], so we utilize these results. The singular part of the free energy approaching 



the critical point depends on p as /sing ~ "^'"S so the operator perturbing the critical theory in 
the integrable direction has scaling dimension Xi^t = 2 — Then we find 

+ 1 in regime 1, 



r+2 
r-1 
r+2 



3 
+ 

in regime 4. 



^ -1-1 in regime 3, 

3 



in regime 2, ^^^^ 



We have written these results in terms of r instead of 9 to emphasize that the derivation only applies 
to r integer, since this is where (14. 2p can be derived. However, we expect that the results can be 
continued to all r within a given regime, since the equations themselves depend on r as a continuous 
parameter. 

A good check of the validity of ()4.3p for generic 9 is that it corresponds exactly to Xi^t = 2s, 
where s is the spin of the discretely holomorphic parafermion ipsiz) described in Appendix B: 
see (iRHjl . 



4.3 Description by conformal field theory 

Here we give formulae for the central charges in all four regimes that are presumably exact. All 
are related to those occurring in completely packed models. However, the doubled loop models are 
not identical to completely packed models: we have checked that the doubled loops have non-trivial 
fractal dimension in regime 4 (see Section [5. 3p . 

One type of critical behavior possible for a doubled loop model is simply to have the two 
colors decouple in the scaling limit. This occurs in the isotropic completely packed version when 
X = —wi = —11)2, and is argued to persist in a region around this point [19]. The central charge 
is simply twice Co(n) given in (13. 6p . Examining our numerical results for c, we see that in regime 
2, the central charge indeed is converging nicely to 2co(n) with the appropriate dependence on n. 
When r is an integer, this is twice the central charge of the conformal minimal models, and so the 
corresponding height model should scale to two decoupled minimal models. An additional check on 
this comes from the fact that the dimension of the integrable perturbing operator in (|4.3p is twice 
that of an operator in a minimal model. Namely, we have X-mt = 2Xi^2, where Xi^2 = 2{r+2) 
scaling dimension of the ^1^2 operator in the minimal model with central charge c = 1 — TfrjYnrr^- 
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It is thus natural to conjecture that in regime 2, the scaling limit of our integrable loop model is 
indeed that of two decoupled completely packed loop models. These conformal field theories have 
been extensively studied [32]. 

In regime 1, the numerics for the central charge are apparently converging to 2cq(„) + ^. Thus 
here the loops apparently decouple as well, but additional critical Ising degrees of freedom appear. 
This is consistent with the mapping to the Ising model at 6 = ^ described in Appendix C. The 
scaling dimension X[^t here is 2Xi^2 + 1; leading to the natural interpretation that the operator is 
a product of the ^1^2 operators in the two minimal models with the energy operator in the Ising 
model, the latter having dimension 1. These extra Ising degrees of freedom, which also occur in 
a certain regime of the square-lattice 0(n) model [33], appear through the following mechanism. 
The vertices of the loop model obey a Z2 symmetry, in the sense that any vertex is surrounded 
by an even number of empty edges. Thus empty edges form polygons where each node has even 
degree, and so they respect the geometry of Ising domain walls for Ising variables a lying on the 
dual lattice. Depending on the values of the Boltzmann weights, these Ising variables may become 
critical in the continuum limit. This is evidently what happens in regime 1. 

The critical behavior in regimes 3 and 4 is not that of two decoupled models. As mentioned 
above in Sec. 13. 1^ in the completely packed version of the doubled loop model, there occurs a coupled 
critical point corresponding to the SUr.(2) x SUr.(2)/SU2r(2) WZW coset model, with central charge 
Cr (I3.10p . The numerical analysis in Fig. [7] nicely fits to Cr in regime 4, and agrees with the Ising 
value c = ^ at 9 = derived in Appendix C. Moreover, when r is an integer, the exponent X-mt 
in ()4.3p belongs to the above coset theory. Thus it is natural to conjecture that the central charge 
throughout regime 4 is c^. As we see from our numerics at n = (see Sec. l5.3p . the fractal dimension 
of a single loop is < 2, so regime 4 represents a "dilute branch" of the coset theory. 

Likewise, in regime 3 the data seem to be converging to + ^, agreeing with the six- vertex 
value c = 1 at ^ = The exponent X^^t is 1 greater than the value in regime 4, so it is natural to 
interpret that the operator is multiplied by the Ising energy operator of dimension of 1. Thus like 
in regime 1, the critical theory presumably includes an extra Ising piece. 

Outside r integer, the conformal field theory in regimes 3 and 4 is not understood. Moreover, we 
will see in the subsequent section that even though the formula for the central charge is applicable 
for all r, it is not even clear whether dimensions of exponents can be continued to values of |n| < 1. 

5 Critical behaviour in regime 4 

The main motivation of this paper is to explore a doubled loop model arising in the truncation of 
the Chalker-Coddington network model. For a connection to disordered systems, the weight per 
loop n and the central charge c must be zero. We have two n = points, but for 9 = ^ inside 
regime 1, the corresponding critical field theory seems to have nothing to do with the CC model. 
Not only do the different colors of loop decouple, but the extra Ising degree of freedom makes c ^ 0. 
We thus in this section focus on the behavior in regime 4, which contains the other n = point at 
9 = ^. 

As noted above, we do not have a conformal field theory description valid in regime 4 outside of 
integer r. It therefore seems a good idea to exploit the fact that the associated height description 
at these points is described by the coset conformal field theory SU(2)r x SU(2)r/SU(2)2r. 

5.1 Integrable perturbations 

This coset theory is known to have two integrable perturbations. One of them, found by using 
level-rank duality on the results of [34], is by the operator with dimension Xi^t = discussed 
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above. This perturbation describes the scaUng Hmit of the height model with elhptic Boltzmann 
weights [25]. In terms of the loop model, we have found that the discrete parafermion ipsiz) of 
Appendix B is the chiral part of the corresponding operator. This parafermion consists of the 
insertion of a one-leg defect for each loop flavour. 

The other integrable perturbation also has a very natural meaning in terms of loops. This 
perturbing operator corresponds to the (l,l;adjoint) operator of dimension ^(i,i;adj) = ^Tpi- Several 
arguments imply that this perturbation corresponds to changing the weight per unit length of the 
loops [33. This integrable field theory describes the scaling limit of an integrable height model 
[36], and using the BWM algebra, it is described in [35] how to relate this height model to a 
dilute doubled loop model very similar to the one we study here. Moving away from the critical 
point in this similar model turns out to be effectively changing the weight per unit length. The 
second argument implying this result involves the S matrices for this integrable field theory, which 
decompose into the tensor product of S matrices of two minimal models Sr x Sr [37] . It is natural 
to interpret the worldlines of a particle in a single minimal model as a loop in the 0(n) model [38] . 
Thus when the S matrix is given by this tensor product, it is natural to interpret the worldlines of 
such particles as doubled loops; when two particles scatter they obey one of the four processes in 
the vertices wi, W2 and x pictured in Fig. |5l In such an interpretation, the weight per unit length 
of the loop is related to the mass of the particle. In the field theory, moving along this integrable 
line corresponds precisely to varying the mass of the particle. 

We denote Xt the thermal exponent, defined as the conformal dimension for the first excited 
state in the zero-leg sector. The numerical calculation of Xt (see Fig. [8]) brings two observations. 
In the region ^ < n, the thermal exponent Xt converges to -^(i,i;adj) = ^rpj even for generic 
values of 9, whereas ^(i,i;adj) was derived only for integer values of r. This indicates that the results 
from the SU(2)r x SU(2)r/SU(2)2r coset WZW model may be continued to arbitrary ^ < 6* < vr. 
However, in the region ^ < 6 < Xt clearly deviates from the continued value this shows 
that not all exponents of the loop model are given by analytic continuation of the WZW coset model 
in this region, including our point of interest 9 = 

5.2 Correlation length exponent u 

The correlation length exponent is defined as the analog for the loop model of uqc (|2-4p . In the 
CC model, the effect of perturbing the energy level E away from the transition value Ec amounts 
to taking /3,/3' out of the critical line ()2.17p . The analog of this perturbation in the integrable loop 
model is to introduce a staggering of the spectral parameter between the even and odd sublattices, 
with symmetric values around the isotropic spectral parameter ipo (j3.19p : = (1 it \)ipo, in the 
range —1 < A < 1. 

The parameter A acts in a similar way to (/3 — (3c) in the original CC model. At A = —1, the 
only allowed loops are those with minimal length, winding around the vertices of one sublattice 
(say, the even one). At A = 1 loops also have minimal length, but wind around the odd sublattice. 
The critical transition takes place at A = 0, where the two sublattices become equivalent, and loops 
may be very long. In the limit A — t- 0, we expect this perturbation to develop a correlation length, 
scaling as 

e~|A|-^ (5.1) 

Since this staggering does not respect the rapidity lines of the square lattice, it breaks integrability. 

Let us first discuss the point 6 = where the model maps to free fermions and remains solvable 
when the A perturbation is included (see Appendix C). At this point, we get the analytical result 
I' = 2, whereas the energy operator of the free-fermion theory is Xt = 1. In Appendix C, we show 
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Figure 8: Thermal exponent Xt in regime 4. Data points were obtained by transfer-matrix diago- 
nalisation, and the solid line represents the exact result for the integrable perturbation dimension 

V — 2r 

-^{l,l;adj) - 744- 



that the effective theory is a massive Majorana fermion with mass proportional to and not A. 
Thus, at 6 = we have the relation between u and the dimension of the perturbing operator: 

Xt = 2-^. (5.2) 

It is natural to assume that both and Xt are continuous in 9, so the effective mass term should 
still be proportional to A^ outside 9 = and the relation (|5.2p holds all along regime 4. 

For 9 exponent u is only accessible numerically, through finite-size scaling. The correlation- 

length exponent v is obtained by assuming a one-parameter scaling law for the energy gap in the 
presence of the staggered perturbation A. For A ~ 0, we expect the behaviour: 



where F is a scaling function. Since eigenvalues are unchanged under A — )• — A, -F must be an even 
function. In particular, at 9 = like for the truncated CC model (see Sec. I2.4p . we get the best 
data collapse (see Fig. [9|) for the value 

z^~l.l. (5.4) 

To our numerical precision, this value may be related through (j5.2p to the thermal exponent Xt — 
1.71. This is an indication that the relation (|5.2p should hold throughout regime 4. 

5.3 Other exponents 

A finite-size scaling plot of the gap corresponding to £i = £2 = 1 is shown in Fig. [TUl The fact that 
it scales to zero faster than 1/L indicates that Xi^i is consistent with zero. This agrees with the 
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Figure 9: Data collapse for the effective thermal exponent Xt{L, A) in the presence of a A pertur- 
bation. The value used for this plot is 1/1/ = 0.9. 



value X|Q|2 = expected for the untruncated model in ()2.2ip . but as far as we know there is no 
fundamental reason for them to agree, and it would be interesting to investigate this further. The 
observed slope in Fig. [TOl suggests existence of an irrelevant operator with scaling dimension ~ 3.2. 

Moreover, we observe numerically that Xt = X2^2- Hence, like for usual dilute polymers, this 
means that Xt is associated to a perturbation of the monomer fugacity (but different from the 
coset-model continuation to r = 0, which would yield -^(i,i;adj) = 0)) and the fractal dimension of a 
path is 

d/ = 2 - ^2,2 ^ 1.71 . (5.5) 



6 Discussion 

Regimes 3 and 4 of the integrable model are particularly interesting, as the two loop colours re- 
main coupled in the continuum limit. They are described by the "dilute branch" of the SU(2)r x 
SU(2)^/SU(2)2r WZW coset CFT (the same CFT as for the completely packed case [H]), with 
an additional Ising degree of freedom in regime 3. Strictly speaking, this theory is only valid at 
the RSOS points, but some critical properties (including the central charge) extend to the loop 
model for generic fugacity n. However, differences between the loop and RSOS spectra exist, as 
shown by our results on the thermal exponent Xt . An analytic study of the loop model through the 
Bethe Ansatz Equations is considered for future work. One also needs to understand if a Coulomb- 
Gas construction (most probably with a three-dimensional target space) could reproduce the coset 
results for generic n. 

The original motivation for the present work was to propose an exactly solvable approximation 
to the IQHE transition. Unfortunately, our results for the correlation-length exponent v clearly 
indicate that the point n = of the integrable model is not in the universality class of IQHE. 
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Figure 10: Log-log plot of the energy gap for the watermelon (1, 1) sector. Data points are fitted 
by a line of slope ~ —2.2, and hence the conformal dimension is Xi^i = 0. 



However, this model is a critical, integrable point in the phase diagram of our modified CC model. 
It should really be considered as the first order in a hierarchy of truncated models, converging to 
the IQHE universality class. Higher-order truncated models should also contain integrable points, 
which may be built by "fusing" the edges of the first-order model, following [36j or a more recent 
approach based on additional Zat symmetry |39j . 

From the numerical point of view, we used the only known efficient method to study a generic 
loop model: transfer-matrix diagonalisation. However, the inherent limitations on the system size 
prevent us from obtaining sharp estimates for the exponents, especially for z/. Recently, new Monte- 
Carlo algorithms have been proposed to simulate 2d loop models [10]. We hope to adapt this new 
approach to two-colour models, and get more precise estimates for the exponents of our truncated 
CC model. 

Acknowledgments The work of P.F. was supported in part by NSF grants DMR/MPS 1006549 
and 0704666. 

Appendix A: The BWM algebra 

In this Appendix, we recall the motivation and definition of the BWM algebra in a graphical 
language. 

The BWM algebra [20] is a braid-monoid algebra, an object relevant to knot theory. It was 
originally designed to compute a certain link invariant, and later it was realised that it could be 
represented by RSOS models related to affine Lie algebras [28]. In [25], a dilute version of the BWM 
algebra was constructed, together with the corresponding i2-matrix. 

In the context of knot theory, the basic objects under consideration are braids. Let (pi, . . . ,pl) 
be L distinct points in the complex plane, and define two copies of each point in three-dimensional 
space, p'j = pj X {0},Pj = pj x {1} , so that the points {p'j} and {pj} lie in two parallel planes. For all 
j = 1, . . . , L, take a curve Tj enclosed between the two planes, and connecting p'^ to p'j. Furthermore, 
impose that the Fj's do not intersect each other. Denote the multiplet F = (Fi, . . . , F^): a braid f3 
is then an equivalence class of F's, modulo continuous deformations of the curves. A typical braid 
is depicted in Fig. [TTl 
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Figure 11: An element of the braid group (left) and a word of the braid- monoid algebra (right) for 
L = 4. 



Multiplication of two braids /3, /?' is defined by the concatenation of the two corresponding 
diagrams, with the convention that diagrams act from bottom to top: the product /5/5' corresponds 
to P above /3'. The /3's form the braid group, generated by the elementary braids Bj,B~^, which 
satisfy the relations: 



BjBj+iBj 
BjBg 



Bn + \BiB 



BgBj 



if li 



> 2. 



(A.l) 
(A.2) 
(A.3) 



Consider now multiplets T = (Fi, . . . ,Tl) of non- intersecting curves connecting all the elements 
of {p'j} U {Pj}, without the restriction that a curve should go from a p'^ to a p'^. The corresponding 
diagrams are then words on the alphabet {Bj, BJ^ ,Ej, j = 1, . . . , (L — 1)}, where the meaning 
of the letters Bj, B~^, Ej is given in Fig. [T2j The algebra on these words is called a braid- monoid 
algebra. It has two parameters {N,oj), and is defined by the braid-group relations ()A.ip - (|A.2p . 
together with the additional relations (see Fig. [T3]) : 



EjEj±iEj 
BjEj 

BiBn + ^E^ 



NE, 
BjBj 



IjJ En 



■'3^j±iJ^j — Ej±iBjBj±i — rjj±iiijj . 
Equations (|A.4p - (|A.5l) mean that the Ej form a Temperley-Lieb algebra with loop weight N . 



E,:+'\ E i 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
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Figure 12: Generators of a braid-monoid algebra. 



The BWM algebra is a braid-monoid algebra ()A.1|) - ()A.7|) where one imposes a linear relation 
between Bj, BJ^ and Ej-. 

Bj — B- ^ uj — 
Ej = l + — ^ , where N = 1 + — . (A. 8) 
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Figure 13: Algebraic rules in the braid-monoid algebra with parameters {N,u}) 



The reason for introducing such a constraint is that the resulting algebra supports a linear form 
(the Markov trace) which is identical to a geometric invariant of the diagrams T (the Kauffman 
polynomial) pO] . 

The dBWM algebra |25] is obtained by allowing vacancies, or equivalently by taking multiplets 
of curves T = (Fi, . . . , F^) with < £ < L. This amounts to adding the generators 



J' 



whose action is depicted in Fig. [TH In equations (lA.ip and (|A.8p . 1 is replaced by Ij, so that the 
Bj, BJ^ , Ej still form a BWM algebra on the set of occupied sites. Additional relations for the 
dilute generators should be included, to implement invariance under continuous deformation of the 
curves in the presence of vacancies. The full set of dBWM relations is given in |25j . 
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Figure 14: The generators of the dilute BWM algebra. 
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Appendix B: Discretely holomorphic parafermion in the loop model 



In this Appendix, we show that the two-colour loop model admits a discretely holomorphic para- 
fermion ipsiz) [411 H2t H3] exactly on the integrable manifold (j3.12p . The parafermion ips{^) is 
defined on the midpoints of the dual lattice C* , and inserts a one-leg defect for each color at point 
z. In the two-point function {^s{0)tps{z)) , there is a black (resp. grey) path 71 (resp. 72) connecting 
and z, and one includes a phase factor involving the winding angles W of the paths 71,72: 

{UO)Mz)) = ^ E E n(C)eT[^(7i)+^(7.)], (B.l) 

(71.72)0^2 C|(7i,72) 

where the first sum is over all possible pairs of paths from to z, the second sum is over the 
loop configurations C compatible with 71 and 72, and n(C) is the Boltzmann weight for a loop 
configuration C. 













/ 


















1 





















Figure 15: A loop configuration contributing to (V's(0)'(/'s(-z))- 



We impose discrete Cauchy-Riemann (CR) on {^ps{0)'^psiz))■ 

Y,{MO)M^)) Sz = , (B.2) 

where the sum is over the edges of an elementary plaquette of C* , and the 5z^s are the corresponding 
elementary displacements. For the discrete CR equations ()B.2p to hold, it is sufficient to fix the 
external loop configuration outside a given plaquette, and ask the total contribution of internal 
configurations to vanish [41 1 H3] . This determines a linear system of equations for the Boltzmann 
weights. To get anisotropic solutions, we consider the analog problem on a rhombic lattice of angle 
a 021 is]. Setting 

X:=e'¥ , ^ := e^°(i+'^) , (B.3) 
we get the 7x7 linear system for the unknowns {t,ui,U2,v,wi,W2,x): 

t + flX~'^Ul — flU2 — V = 

ri^ui — \~'^U2 — iJiX^v + ^\~'^{v?wi+W2 + 2nx) = 

— A^tii + n'^U2 + fiX'^^v — fi{wi + v?W2 + 2nx) = 

-^iX^ui + /iA-^n2 + n'^v - X^wi - X-'^W2 - 2x = (B.4) 

nXui — nX~^U2 — fJ.X~^v + fj.X~^[n{wi + W2) + {'n? + l)x] = 

7T,A~"^Mi — An2 — n^Xv + ^X~^{nwi + x) + fJ,X{w2 + nx) = 

—ui + nX'^U2 + nfiX~'^v — iJ.X~'^{wi + nx) — fiX'^{nw2 + x) = 0, 



26 



As a first step, we need determine the spin s by going back to the isotropic case a = tt/2. Imposing 
ui = U2 and wi = W2, ()B.4p reduces to a 5 x 5 system, whose determinant is: 

-X-\X^ + l)(nA^ - 1)\X^ + A-^ + - 3n) . 



D{n, A) 

Using the parameterisation n = —2 cos 29, this determinant vanishes when: 

exp(2i7rs) = exp{±6i9) . 
For a general angle a, the solution of ()B.4p is a set of a-dependent weights t{a) 



apply the substitution: 



a 



1 + s 



(B.5) 

(B.6) 

, x{a). If we 

(B.7) 



we observe that the solution of (jB.4p is identical to the integrable weights (|3.12|) . This is analogous 
to what was found for various other integrable models with a discrete holomorphic parafermion |4H 
321 133]. Note that the relation ()B.7p is consistent with the discussion on the Fermi velocity in 
Section 13.31 

Moreover, at the Ising points ^ = f i x (^^^ Appendix C), the spin s = ^ is consistent with (|B.6p . 



J. 

(1) 



r \ 



(2) (3) (4) (5) (6) (7) 

Figure 16: External loop connectivities outside an elementary plaquette. 



Appendix C: Free fermions at ^ = |, ^ 
C.l Mapping to a staggered 8V model 

At ^ = ^ and = the model maps to a free fermion Hamiltonian through the mapping sequence: 

two-color loop —7- square-lattice 0(n = 1) — )• staggered 8V — t- free fermions. 

In both cases, we start from the two-color loop model and do the sign change (n, x) — t- (— n, —x) 
to get n = 1. With this value of ?i, one may discard loop connectivities, and the model is simply 
local, with occupied and empty edges. After the exchange of occupied/empty edges, we get a 
square-lattice 0(n) model |33j with n = 1 and weights : 

t = Wl+W2 — 2x, Ul = Ui, U2 = U2, V = V, Wl + W2 = t . (CI) 

(The tilde is here to avoid confusion between the 0(n) and two-colour loop model Boltzmann 
weights). For = we have the specific values 



t 


= sin 2ip — "v/S 


Ul 


= -v/S cos if 


U2 


= — \/3 simp 


V 


= sin 2lp 


Wl 


= — cos(^ — 2ip) 


W2 


cos(f + 2lp) 



(C.2) 
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which are exactly those of the integrable square-lattice 0(n) model [33] with n = 1 at the dilute 
critical point. The 0{n = 1) model maps in turn to an eight-vertex model (see Fig. [T7|) with 
staggered weights 



UJ3 



UJ2 



(W5,W6) 
UJ'j = 



\/3 cos 93 
— \/3sin(^ 

{(— sin 293 — \/3, sin 293 — 1/3) on even sites 
( sin 2lp — -v/S, — sin 2(^9 — \/3) on odd sites 
sin 2ip . 



At 



one gets the same 8V model, up to irrelevant signs. 



UJ2 



W3 



UJ4 



UJ5 



< > 



W7 



+ + 



(C.3) 



COS 



Figure 17: Correspondence between the 8V model and the 0(n = 1) model. The mapping depicted 
here is valid on one sublattice, say even sites. On odd sites, all arrows must be reversed. 



C.2 Very anisotropic limit: the XY chain in a magnetic field 

In terms of the Pauli matrices aj, the 8V ^-matrix reads: 

+ {uJ3 (^Jcr-^+i + W4 0-/0-7+1 + W7 crJcrJ+1 + ^8 o-+f7+.i) 

- uj2){a^ + (t|+i) + ^(ws - u;6)(ct| - . (C.4) 

We can now take the very anisotropic limit 93 — )• 0. Denoting by a prime the derivative with respect 
to 93 at (/3 = 0, and using the weights ()C.3p . we obtain: 

i?;8V ^ -Vsiajaf^, + alaj^,) + 2iaTaJ^, + - - . (C.5) 

The critical Hamiltonian is given by 

where Jx = {2 — \/3), = (2 + \/3) and h = 2. The alternating sign of the last term in ()C.6P can 
be eliminated by the unitary change of basis: 

L/2 

Hq^U^HoU, where U :=Yla^f,_^ . (C.7) 

^=1 
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This maps Hq to an XY chain in a magnetic field [44^ 145] : 

L 



H^Y = U^HoU = - i ^ [(1 + ^) a]a]^, + (1 - 7) + 



(C.8) 



where 7 = — ^ and h = 2. In particular, we have learnt that this particular point of the XY spin 
chain is exactly equivalent to the integrable dilute 0(n = 1) model. 

C.3 Staggered perturbation associated to u 

The setting of the 8V model also allows us to consider a staggered perturbation like the one defining 
exponent (see Sec. 15. 2p . To do this, we introduce staggered spectral parameters (1 ± \)ip on the 
even/odd sites. We take the very anisotropic limit (/3 — )• 0, with A fixed. This way, the parameter A 
controls the strength of the perturbation. The resulting Hamiltonian has the form 

Hi\) = Ho + XHi 

where 

j=i i=i 

and Jx,Jy are the same as for Hq. After the unitary change of basis defined by U, we get the 
perturbing term: 



A 



H, = U^H.U = -- J2i-iy (1 + 7) + (1 - 7) cT]al 



+1 



(C.IO) 



C.4 Exact free-fermion solution 



In this paragraph, we expose the exact solution of H{\) for general values of 7, h, A. Following 
we can solve the model H{X) by a Jordan- Wigner transformation, mapping the Pauli matrices aj 
to fermion operators 



0-1 

n 



0-1 

n 



CTf (7- 



J, 



1 



c;c. = -(l-a|; 



"3-3 2 



(C.ll) 



obeying anti-commutation relations: 

{cj,ce} = , {cj, 4} = 6je . 
In this language, the perturbed Hamiltonian reads 



(C.12) 



j=l { 



t J 



We introduce two species of fermions 



C21 '■— C21-1 



(C.13) 



(C.14) 
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and their Fourier modes 



^ L/2 
V £=1 



27rm L 
— ; TT, m = l,..., — . 

L/2 ' ' ' 2 



(C.15) 



We can now rewrite H{X) as 

— 7r<(jr<7r M,!' 



(C.16) 



where the matrices Aq and Bq read 



^ a* 





/3, 



:= -[(1 + A)e^'' + (1-A)], 
fiq := 7[(1 + A)e*'?-(1-A)]. 



(C.17) 



Like in [H], the energies e^{q) are the square roots of the eigenvalues of {Ag + Bg){Ag — Bq): 



ei,2(g) = 2l/^ + (1 + 72 A2) cos2 I + + A2) sin^ | ^ ^/472A2 + 



cos — h iA sm — 
2 2 



. (C.18) 



This is the two-branch dispersion relation for arbitrary 7, /i, A. 

To express the corresponding eigenmodes, we need the unitary 2x2 matrices Wg, Vg defined by 
the linear relations 



(Ag - Bg)Wg = Vg D g , (Ag + B g)Vg = WgDg , Dg 1 = 

The Bogoliubov transformation diagonalising H{X) is 



e2{q) 



^E[(^ + ^) 



tJ.u,q (^u,q 



{W -V) 



^iu,q ^u, — q 



Unitarity of Vg and Wg ensures the canonical anticommutation relations 
In terms of the r/'s, the Hamiltonian reads 

— 7r<g<7r /j=l,2 



(C.19) 



(C.20) 



(C.21) 



(C.22) 



Note that there are L/2 distinct momenta q, and that each momentum corresponds to two modes 
^ = 1,2. Thus, we recover L independent modes rj^^g. 
As a final step, we perform the change 



if ^ > 0, 
'ni,q if g < 0, 



so that the modes with (7 < are now considered as holes, the Hamiltonian becomes 
HW= E ^^M)^^^,q^h,q^ where := sgn(g) e^(g) . 

— 7r<g<7r /^=1,2 



(C.23) 



(C.24) 
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C.5 Critical Majorana fermion at h = 2, X = 

The dispersion relation of the XY chain in a magnetic field is obtained by setting A = in (jC.lSp : 



61,2(9) = 2y cos +72sin2| 

For h = 2, ei is critical at g = 0, whereas €2 is not criticaQ: 



(A = 0) 



(C.25) 



ei g = 4sin-4/sin — h 7 cos^ - 
4V 4 4 



£2(9) = 4sgn(g) cos cos^ ^ + 7^ sin^ ^ 



(C.26) 
(C.27) 



The dispersion relation ^i{q) is approximately linear at q = (see Fig. llSp . In the ground state, all 
levels with — vr < q < are filled: this is a Fermi sea with only one Fermi level qf = 0, and thus it 
corresponds to a Majorana fermion with central charge c = ^. The critical eigenmodes are obtained 
from (lU^Ol) : 



^(ci,, + C2,,) + isin%^(4^_^ + 4 



V2,q 



where 



'ni,q — COS — v<-l,q -r ^2,97 -r '^iii — V-l.-g "T , 

COS %^ (ci,g - C2,q) + i sin %^ (cj - 4 _^) 



:= Arg (1 — cos I — i sin I) , 



Arg (1 + cos f + i sin |^ 



(C.28) 



(C.29) 



After the change r] ^ r] (IC.23p . the continuum limit is described by the effective Hamiltonian 





Figure 18: The dispersion relation for /i = 2,7 = — Left: critical case A = 0. Right: A 
The full (resp. dotted) lines represent ?i (resp. e2)- 



^ In the case h — j = 0, both modes ei,e2 are critical at g = vr, and the corresponding field theory is a critical 
Dirac fermion. 
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C.6 Gapped theory at h = 2, X > 

When the A perturbation is turned on, an energy gap opens at q = (see Fig. [T8|) : 

AE ■.= 2ei{q = 0) = 2 ^/TTj^ -1 ^ J^X" , v = 2. (C.31) 

To understand this value of v, we shah analyse the perturbing term Hp in terms of the critical 
modes r/i^g (|C.28p . The relations (IC.28P can be inversed, to give: 



(C.32) 

Cl,g - C2,g = COS r]2,q - I siu ^ r]J^ _g := a2,g ■ 



ci,g + C2,g = COS ^ - z sm ^ r/l _g := ai.g 
^2.. ^„ ^2,, ^t' 



The perturbing term has the expression 



^P = l Yl { (1 - cos q) {a\ ga^^g - 4,g"2,q) 



-7r<g<7r 



2 



+1 sm g 

-7(1 + cosg) (al o{ +ai _qa2,g) }. (C.33) 



In the region q 0, the first term in ()C.33P is of order g^, and thus it generates irrelevant terms 
of the form 'rj\d'^rj^ in the continuum limit. The second term is of order q, and corresponds to 
rj^dxTJi, which renormalises the Fermi velocity. At first order in A, the third term has no effect on 
the continuum theory. However, in second-order perturbation in A, it generates terms of the form 
(^2^2)(^i^i)' which are non-vanishing since the lowest rj2 modes are occupied in the ground state. 
From this analysis, we obtain the effective Hamiltonian in the continuum limit 

ff.,(A)oc/d. (i,ia.5. + const x,^A' ■qlrjA . (C.34) 

The "mass term" has dimension Xt = 1, and the energy gap thus scales as 

AE oc (A^)^ . 
Comparing with ()C.3ip . we get the scaling relation 

Xt = 2-^. (C.35) 
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